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ABSTRACT 

Aims. We investigate the spectral and temporal behavior of the high mass X-ray binary Vela X-1 during a phase of high activity, with 
special focus on the observed giant flares and off states. 

Methods. INTEGRAL observed Vela X-1 in a long almost uninterrupted observation for two weeks in 2003 Nov/Dec. The data were 
analyzed with OSA 7.0 and FTOOLS 6.2. We derive the pulse period, light curves, spectra, hardness ratios, and hardness intensity 
diagrams, and study the eclipse. 

Results. In addition to an already high activity level, Vela X-1 exhibited several intense flares, the brightest ones reaching a maximum 
intensity of more than 5 Crab in the 20-40 keV band and several off states where the source was no longer detected by INTEGRAL . 
We determine the pulse period to be 283.5320 + 0.0002 s, which is stable throughout the entire observation. Analyzing the eclipses 
provided an improvement in the ephemeris. Spectral analysis of the flares indicates that there appear to be two types of flares: relatively 
brief flares, which can be extremely intense and show spectral softening, in contrast to high intensity states, which are longer and show 
no softening. 

Conclusions. Both flares and off states are interpreted as being due to a strongly structured wind of the optical companion. When 
Vela X- 1 encounters a cavity with strongly reduced density, the flux will drop triggering the onset of the propeller effect, which inhibits 
further accretion, giving rise to off states. The sudden decrease in the density of the material required to trigger the propeller effect in 
Vela X-1 is of the same order as predicted by theoretical papers about the densities in OB star winds. A similarly structured wind can 
produce giant flares when Vela X-1 encounters a dense blob in the wind. 

Key words. X-rays: stars - stars: flare - stars: pulsars: individual: Vela X-1 - stars: magnetic fields 

1 . Introduction The neutron star revolves wi t h a long spin period of ~283 s 

«i v , n„ n „ m ,m ■ ,• u- u v u- (RaDDaDo rt&"McChntock|[T97l [McClinto ck et al.lll976l) . Both 

Vela X-1 (4U 0900-40) is an eclipsing high mass X-ray bi- ^ spin period and spin period derivative have changed elTat . 

nary (HMXB) consisting of the B0.5Ib super giant HD 77581 ically since their first measurements, as expected for a wind- 

and a neutron stai _xyith an orbital period of 8.964 days accreti ng system. The evolution of the spin pe riod is most ap- 
dvan Kerkwiiketal.1 [1995J) at a distance of ~2.0kpc <|Nagase| priately deS cribed by a random walk model (Tsunemi 1989; 

11989J). The optical companion has a mass of -23 M Q and a radius ziolkowski 1985) . Although the source exhibits strong pulse-to- 
of -30 R Q d vanKerkwiiketaU | 1995|) . Due to the small separa- ke variationSj a pu i se . profi i e fo i ded over several pulse periods 

tion of the binary system with an orbital radius of just 1 .7 R*, the shows remarkable stabilit (for 10 pulses or more; staubertiHD 

massive 1 .9 M neutron star (Vela X- 1 is the most massive com- [T^ eyen oyer decadefj ( Rauben heimi[l9l|. At energies be- 

pact object known to be a neutron star; | Quaintrell et alj | 2003 fc low 5 key ±e pulse . profile consists of a complex five . pe aked 

| BarzivetalJ | 200l|) is deeply embedded in the dense stellar wind structurei which transforms at energies above 20 keV into a sim- 
of its optica companion HD 7758 1 (M* = 4 x 10" 6 M yr 1 ; k double . pe aked pulse-profile ( Stau birtltd] fT980) where the 

| Nagase etalj | 1986|) . X-ray lines indicate that this wind is mho- tWQ kfj m ± ht tQ be due tQ ^ tWQ accreting magnetic 

mogeneous with many dense clumps dOskinova et al.1120071) em- po i es of the neutron star 
bedd ed in a far thinner, highly ionized component (Sa ko et all 

119991) " With ~4 x 10 erg s , the X-ray luminosity of Vela X-1 

is typical of a high mass X-ray binary. Observations in the 

Send offprint requests to: I. Kreykenbohm, past, however, have shown that the source is variable with ob- 

e-mail: Ingo.Kreykenbohm@sternwarte.uni-erlangen.de served flux reductions to less than 10% of its normal value 
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(off states; IKreykenbohm et all [l999t iKretschmar et al.1 119991: 
llnoue et all 1 19841) . while periods of increased activity have also 
been observed during which the flux increases within an hour 
to a multiple of th e previous value, reaching peak flux lev- 
els cl ose to 1 Crab (IKreykenbohm et a i]|l999b lHa berl & White! 
ll990tlKendziorraetalJll989l) . In this respect. Vela X-l is sim- 
ilar to sources such as 4U 1700-377 and 4U 1907+09, for 
which low luminosity states and flares have also been ob- 
serv ed, as is rather typical for wind-accreting systems (see 
e.g. iFritz et al.1 120061: Ivan der Meer etai] 120051: Tin 't Zand et alj 
119971: lHaberl et al.l Il989t ). Although Vela X-l is a well stud- 
ied object, only observations by INTEGRAL revealed that the 
flares in Vela X- 1 can be brighter than p reviously anticipated 
dStaubert et alJl2004tlKrivonos et al.ll2003l) . While the flaring ac- 
tivity is thought to be due to a strongly increased accretion rate, 
M, the origin of the M variations is unknown. 

The phase averaged X-ray spectrum of Vela X-l was usually 
modeled with a po wer law modified at higher ene rgies by an 
exponential cutoff dTanakalll986l:IWhite et al.lll983j or with the 
Negative Positive Exponential (NPEX-model: lMihar3 Il995l) . 
The spectrum was found to be modified by strongly orbital- 
phase-dependent photoelectric absorption at lower energies 
due to the dense s tellar wind and an ac cr etion wake trailing 



the neutron star ( Goldstein et a l] |2004|: Kreykenbohm et al] 



1999; iFeldmeier et alJI 1996b lHaberl & Whitdl 19901) . At 6.4 keV. 
an iron fl uorescence line and occasionally an iron edge at 
7.27 keV dNagase et alj Il986l) were observed in the X-ray 
spectrum. At higher energies, cyclotron reso nant scattering 
featur es (CRSFs) b etwee n 25 and 32 keV ([M akishim a et alj 
1 19921 IChoietalJ 119961; iKrevk enbohm etaij |2002|) and 



at ~55keV dKendziorra et alj 119921; lOrlandini et all 119981; 
Kreykenbohm et al.1 119991: lLa Barbera etai] 120031: lAttie et alj 
20041) were present, although the interpretation o f the 25keV 



feature is still sometimes debated dOrlan dini 2006). 

The remainder of this paper is structured as follows. In 
Sect. [2] the data and software used are described. Section [3] de- 
scribes first the temporal analysis of the data, i.e. light curves, de- 
termination of the pulse period, quasi-periodic modulations, and 
then the analysis of the eclipse and the spectral analysis. The 
results are discussed in Sect. [4] and a summary is presented in 
Sect.g] 



2. Data 



2.1. Instrument and Software 



The INTEGRAL observatory dWinkler et al.l2003l) is in a highly 
eccentric orbit with a period of 71 h 49 m , ideal for long unin- 
terrupted observations with a low X-ray background. Due to 
the high eccentricity of the orbit the perigee passage (when 
INTEGRAL is inside the radiation belts) has a duration of only 
~8 h, which minimizes the time when no science observations 
are possible due to the high radiation background. 

INTEGRAL has four science instruments, which provide 
coverage from 3keV up to lOMeV as well as in the opti- 
cal: the imager IBIS/ISGRI dUbertini et alJ I2003L 20keV to 
800 keV) with moderate energy resolution and a large effective 
area, the spectrometer SPI (20keV to lOMeV; IVedrenne et alj 
|2003[) for the analysi s of nuclear lines, the X-ray monitor JEM- 
X (lLundetai]|2003l 3keV to 35 keV), and the optical moni- 
tor OMClMiSsiSiiEiDllQQl)- Al l high-energy in struments 
are coded-mask telescopes (see e.g. lln'tZandl II 992, for a re- 
view of this technique). To improve the imaging quality, the 
satellite performs raster observations of the vicinity of an X- 
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Fig. 1. Intensity mosaic image using ISGRI data of the Vela re- 
gion. All five revolutions from Rev. 137 to Rev. 141 have been 
used for this mosaic. Vela X-l is by far the most significant 
source. Several more sources are also detected. Note that the 
times of the eclipses of Vela X-l (see Fig. [2j have been excluded 
from this mosaic, resulting in an exposure of ~960ksec. The 
noisy rim of the mosaic is due to the low coding factor in the 
outermost part of the partially coded field-of-view. 



ray source, retaining the target in the field-of-view of ISGRI. 
Due to this "dithering" strategy the off-axis angle of the target 
source changes significantly during the observation. When the 
source is more than ~4°5 away from the pointing direction, it is 
in the partially coded field-of-view of ISGRI. With increasing 
distance from the pointing direction, the coding factor decreases, 
causing increased uncertainties in the images, flux values, and 
spectra. Individual pointings made during these dithering obser- 
vations are called science windows (SCWs). They have typical 
durations of 1800 s, 2200 s, or 3600 s. These SCWs are then asso- 
ciated with INTEGRAL revolutions, i.e. complete orbits of the 
INTEGRAL satellite around the Earth. 

To prepare the INTEGRAL data for analysis, we used the 
Offline Science Analysis Software (OSA), version 7.0, and its 
associated calibration files. In particular, we make extensive 
use of the tool HJight. We carefully checked the behavior of 
HJight (see AppendixlAl, because the IBIS cookboolQ cautions 
that HJight should only be used to analyze the timing behavior 
within a given science window. For further analysis, we used 
HEADAS release 6.2. Spectral fitting was done with XSPEC 
11.3.2ad dDorman & Arnaul200 lb LArnaudll 19961) . 

2.2. Data 



As part of the AOl core program dWinklerll200l . INTEGRAL 
observed the Vela region (see Fig. [TJ continuously for five con- 
secutive INTEGRAL revolutions from the beginning of rev- 
olution 137 (JD 2452970.86) until the end of revolution 141 
(JD 2452970.86) resulting in approximately 1 Msec of data (see 
Fig . [2j» . The observation was performed in a 5 x 5 dithering pat- 
tern with stable pointings 2° apart. 

We chose to use all available science windows from Rev. 137 
to Rev. 141 to be able to derive a contiguous light curve with as 
few interruptions as possible (data gaps due to the perigee pas- 
sage of the satellite are obviously unavoidable; see Fig. [2]). Since 



1 available at 
http : //isdc . unige . ch/?Support+documents 
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Vela X-l is a bright source, the OS A software has no problem 
in detecting the source and determining its flux level accurately 
even when the source is at an off-axis angle of more than 14°; in 
any case, fewer than 5% of the pointings had an off-axis angle 
larger than 14°. For studying the timing behavior on timescales 
short compared with a SCW, i.e. period determination and search 
for QPOs (see Sect. 13.21 ), the absolute flux is not important and 
the temporal properties are unaffected by a non-optimal off-axis 
flux correction. 

Figure Q] shows the image of the Vela region from these ob- 
servations. While Vela X-l is by far the brightest source in the 
field-of-view of ISGRI, we also detect 4U 0836-429 as a very 
prominent source reaching about a third of the average inten- 
sity of Vela X-l, and the relatively weak sources H 0918-5459, 
the Vela Pulsar, and two sources first reported by INTEGRAL 
dden Hartog et alj l2004t ISazonov et ail 120051) . Since the two 
brightest sources Vela X-l and 4U 0836-429 are well separated 
(about 6°7) and all the weaker sources are even more distant, con- 
tamination of the spectrum of Vela X- 1 due to the presence of the 
other sources is of no concern. Data from JEM-X and SPI have 
not been used in this analysis due to the far smaller field-of-view 
of JEM-X and since Vela X-l is off-center in the observed field 
(Vela X-l was only within the fully coded field-of-view of JEM- 
X for less than ten out of the ~550 individual pointings). The 
SPI instrument, on the other hand, provides a high spectral reso- 
lution, although, due to its low effective area, it is not possible to 
study data on timescales of seconds as required here. 

3. Data analysis 

3. 1 . Light curves and flux 

Vela X-l was found in a strongly variable state during the 
Nov/Dec 2003 observation by INTEGRA L While periods of 
increased activity were observed before (iKrev kenboh m et all 
ll999HHaberf 1994), the behavior found in this observation (see 
also Staubert et al. 2004) is indeed extreme. 

Most prominently, on 2003 November 28 (JD 2452971.67), 
INTEGRAL observed the brightest flare ev er seen from 
Vela X - 1 (designated flare 1; see Fig.|2]and also iRrivonos et ail 
120031) . During the flare, the 20^-0 keV count rate increased from 
a SCW averaged pre-flare value of ~55 counts s (~300mCrab, 
or 1.6 x 10~ 9 ergs~') by a factor of more than seven to 
405 counts s (2.3 Crab) within only 90 minutes - normal flar- 
ing a ctivity reaches peak flux values not higher than 1 Crab (see 
e.g. IKrevkenbo hm et all fl999l) . Vela X-l was therefore more 
than ten times brighter than on average (see Table [1]) which im- 
plies that flare 1 is a giant flare. The average flux values for 
Vela X-l in Table Q] were derived from a spectrum with 400ksec 
of exposure obtained between revolution 8 1 and 89 in June/July 
2003, when Vela X-l was in a normal state. In the following, we 
designate a flare with a peak flux of more than 2 Crab as a giant 
flare, as opposed to normal flares, which do not reach this flux 
level. These flares reach their peak rapidly, i.e. r r ise/r toto i < 0.3, 
where r,i se is the time from the onset of the flare to the peak. 
While flare 1 was detected in all energy bands, it was most pro- 
nounced in the 20-30 keV band. As stated above, Vela X-l was 
never in the field-of-view of JEM-X during the flare, such that 
coverage at even lower energies is unavailable. 

After the peak (duration about half an hour, see Fig. [3]), the 
flare decayed quickly in less than 2h to an intensity level of 
< 1 Crab and within ~ 1 1 h to a SCW averaged count rate of 
~35 counts s _1 (200 mCrab), somewhat lower than before the on- 
set of the flare (see Fig. [2J. The total energy released between 



Table 1. Flux values of Vela X-l in science window 
013700420010, during which giant flare 1 reached its maximum. 
The increase indicates the factor by which Vela X-l was brighter 
during the flare than during its normal state. The flux values are 
averages for the entire science window. For the peak flux reached 
in the pulses, see Table [2] 



energy 


Flux (normal) 


Flux (flare) 


Flux (flare) 


increase 


[keV] 


[ 10-'° erg 


s cnT 2 s _1 ] 


[Crab] 




20 - 30 


10.7 


136.3 


2.8 


13 


30 - 40 


5.6 


58.6 


1.9 


10 


40 - 50 


1.9 


16.6 


0.7 


9 


50 - 60 


0.6 


4.5 


0.3 


8 


60 - 80 


0.7 


2.6 


0.1 


4 


80 - 100 


0.1 


0.6 


<0.1 


6 



Time [MJD-52971] 
0.14 0.15 0.16 0.17 0.18 0.19 




2003 November 28 



Fig. 3. Close up of the light curve with a time resolution of 20 s 
of the peak of giant flare 1 on 2003 November 28. The peak is 
reached at MJD 5297 1.18 with 923 counts s~ 1 (corresponding to 
5.2 Crab) in the 20 to 40keV band. This value is significantly 
higher than the count rate averaged over one science window, 
since the count rate is not constant. Furthermore, the source is 
strongly pulsating, exhibiting the well known double pulse. 



20keV and 40keV amounts to 1.15 x 10 41 ergs. We empha- 
size that these fluxes and those given in Fig. [2] and Table Q] are 
SCW averaged fluxes. These values average over fluctuations on 
shorter timescales such as pulsations, during which the source 
sometimes reached significantly higher intensities (see below 
and Fig. [3]). 

In the following days, three more flares (flares 2 to 4, see 
Fig.|2|, as listed in Tabled were observed. All three flares were 
shorter and less intense than flare 1 on a science window av- 
eraged basis, but still reached SCW averaged intensities close 
to 1 Crab. Unfortunately, INTEGRAL was in engineering mode 
during flare 3 and the standard OSA pipeline rejects these data; 
therefore, only count rates could be obtained and no further anal- 
ysis of this flare is possible. 

On 2003 December 7 (JD 245981.10), another intense flare 
was observed (designated Flare 5, see Fig. |2). Unlike flare 1, 
during which the brightness of the source increased rapidly, it 
took ~8 h for flare 5 to reach its SCW averaged maximum 20- 
40keV flux of -1.2 Crab. The decay lasted ~5h until Vela X-l 
reached its pre-flare count rate of ~35 counts s _1 (200mCrab in 
20-40 keV). Although quite bright, flare 5 is therefore signifi- 
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Table 2. Overview of the observed flares. See Fig[2]for the num- 
bering of the flares. The time is the onset of the flare. To obtain 
the peak fluxes F pea k, a light curve with a time resolution of 20 s 
was used. 7ri se is the time from the onset of the flare to the peak, 
while Ttotai 1S the duration of the flare. 



Flare 


Time 
[MJD] 


Duration 
M 


[Crab] 


r ri , c 


Remarks 


1 


52971.15 


11200 


5.2 


0.15 


giant flare, spectral 












softening 


2 


52975.34 


5 200 


2.6 


0.83 


no spectral change 


3 


52976.50 


1800 


5.3 


0.28 


giant flare, very short 


4 


52977.15 


12 900 


1.9 


0.13 


spectral softening 


5 


52980.31 


31400 


3.9 


0.63 


high intensity state, 












no spectral change 



Time [MJD-52976] 
0.54 0.55 0.56 0.57 0.58 




2003 December 3 



Fig. 4. ISGRI 20^10 keV light curve of the short flare 3 on 
2003 December 3 with a time resolution of 20 s. The count 
rate increases from <40 counts s at the onset of the flare to 
929 counts s _1 during the peak (thus exceeding even flare 1) in 
less than 800 s only to decay within 1200 s to a count rate of 
<60 counts s _1 again. Since the flare is extremely short, its peak 
intensity is not evident in a light curve using SCW long bins as 
in Fig. |2 

cantly less intense than giant flare 1, and also far longer, i.e. it 
is a high intensity state. In the following, we designate flares 
that have a duration of many hours and exhibit a rather slow in- 
crease, i.e. r r i se /r to t a i > 0.5, as high intensity states, as opposed 
to flares which exhibit a rapid increase, such as flare 1. With 
1.22 x 10 41 ergs the fluence of flare 5, however, was very similar 
to the energy release of flare 1 . 

A more detailed analysis using a light curve with a 20 s time 
resolution showed that the source pulsated all the time, including 
during flares, and exhibited the well known strong pulse-to-pulse 
variations (Pig.[3j. In the maximum, a peak 20^-0 keV count rate 
of ~920 counts s _1 in one 20 s time bin (5.2 Crab) was observed 
for flare 1 and ~930 counts s _1 (5.3 Crab) for flare 3. Flare 3 
on December 3 was therefore also a giant flare (see Fig. |4|. 
However, flare 3 was significantly shorter: the entire flare lasted 
less than 2000 s, but it was as bright as flare 1 (see Table|2]l. 

Extending the analysis to the non-flaring parts of the light 
curve, we detected a quasi-periodic oscillation (QPO), similar to 
other accreting X-ray pulsars (e.g. in 4U 01 15+63 or V 0332+53 ; 
iHeindl etalj[l999i: iMowlavi et al]|2006t iTakeshima etaT]| 19941) . 



Table 3. Observed off states during our Vela X-l observation. 
The duration is the time between the onset of the off state (i.e. 
sudden intensity decrease) and end of the off state (sudden in- 
tensity increase; see Fig. [6|. Off state 5 is special as the source 
stays at a very low luminosity level after the end of the off state; 
in this case, the duration equals the time when no pulsations are 
detected (cf. Sect.|4~3ll. 



Off state 


start 


stop 


duration 




[MJD] 


[MJD] 


[s] 


1 


52981.095 


52981.105 


860 


2 


52981.275 


52981.291 


1380 


3 


52981.348 


52981.365 


1470 


4 


52981.373 


52981.379 


520 


5 


52981.422 


52981.445 


1980 



Detecting a QPO, especially if it is transient, can be difficult. 
Using dynamic power spectra (PSDs), we found a short-lived 
QPO with a period in the range of a few 1000 s, but no evidence 
for the presence of any quasi-periodic behavior below 140 s (half 
the pulse period, Sect. 13.21 1. The short-lived QPO with a period 
of ~6820sec appears to be regular and inconsistent with pure 
stochastic behavior (see Fig. [5}. Subsequent period searches on 
the corresponding data subset clearly detect the period. We note 
that the quasi-periodic modulation shown in Fig.|5]is far stronger 
and inconsistent with the NOMEX effect, which can cause in- 
tensity variations from SCW to SCW, but not within a given 
SCW. 

In this case, the amplitude of the oscillations was large, i.e. 
the mean count rate was -45 cps with an amplitude of 25 cps. 
During the bottom part of the oscillation, however, the flux some- 
times decreased to approximately zero counts for a short time 
(see Fig. |5), i.e. the source turns off (see the following discus- 
sion). 

Although these modulations are evident (see Fig. |5), the 
statistics for the entire observation are weak since this event was 
short lived. Furthermore, it is well known that pure red noise can 
also produce quasi-periodic flu x variations (for a dis cussion us- 
ing Mrk 766 as an example, see iBenlloch et al. I I200T1) . therefore 
these modulations must be treated with care. 

During and after the part of the light curve in which the 
quasi-periodic oscillation was present, we observed several off 
states, during which no significant residual flux was detectable 
by ISGRI and no modulation due to the otherwise omnipresent 
pulsations was present (Fig. |6}. The onset of these off states oc- 
cured without any identifiable transition phase. The source sim- 
ply switches off, or more precisely, the luminosity of the source 
drops below the detection limit of ISGRI. The same statement is 
also true for the end of the off states, where Vela X-l switches 
on again and immediately resumes its normal (pre off state) in- 
tensity level and behavior. We identified five such off states dur- 
ing the entire observation, which all occured in the 12 h from 
MJD 52981.0 and MJD 52981.5 (Table 0. Figure shows off 
states 3 and 4, which are separated by only 600 s (about two and 
a half pulse periods). In between the two off states, Vela X-l pul- 
sated normally and at a normal intensity level. Since the two off 
states occured in rapid succession, they could also be considered 
as one longer off state interrupted by 2 normal pulses. 

Off state 5 differed from the other four off states: the onset 
of off state 5 is also sudden, but after the pulsations are observed 
again, the flux level remains low for several thousand seconds. 

2 For an expl anation of the N OMEX sup port structure and i ts impli- 
cations, see e.g. lLubinski et all i2004h and[Reglero et aT] d200lh . 
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Time [MJD-52980] 
1.1 1.2 



c 

3 

O 

B 

H 
C 
3 

O 

U 




23 01234567 
2003 December 8 



Fig. 5. A closeup of a part of the pulse averaged light curve of 
Vela X-l (i.e. a light curve with a time resolution of 283 s). A 
quasi-periodic modulation with a period of ~6820s is evident 
in this part of the light curve, as indicated by the overplotted 
sine wave. Note that during the trough between 2 h and 3 h, and 
especially following the quasi-periodic modulation, the count 
rate decreased several times to zero for a short time, i.e. the 
source became completely undetectable by ISGRI (see Table [3] 
and Sect. l3.1l for discussion). 



Time [MJD-52981] 
0.35 0.36 0.37 0.38 



0.34 




8:30 8:40 8:50 9 

2003 December 8 

Fig. 6. 20^40 keV light curve with 20 s resolution of off states 3 
and 4 of Vela X-l, which are most remarkable as the source sud- 
denly becomes undetectable by ISGRI and then turns on again 
two times within one hour. During the off state no pulsations are 
discernible, while outside the off state the double pulse-profile 
of Vela X-l is visible. 



Off states 1 to 4 are shown in Fig. [5] where the flux declines to 
almost zero for a short time (note, however, that the light curve 
in Fig.[5]has a different binning). 

3.2. Pulse Period 

Due to strong pulse-to-pulse variations, the determination of the 
period of pulsars with long pulse periods such as Vela X-l is 
rather difficult. We th erefore used a two step approach. First, 
using epoch folding (Leahy et all 1 19831) . we derived an ap- 



proximate period of 283.5 ± 0.1s for the Vela X-l observa- 
tion. We then used this period as a starting point for a pulse- 
profile template-fitting approach. In this approach we derive 
sets of pulse-profiles from the beginning to the end of the ob- 
servation using sufficient data for the long-term pulse-profile 
to emerge. Since the lon g-term pulse-profile i s known to be 
extremely stable (see e.g. IStaubert et alJ Il980t iRau benheimer 
1990; iRrevkenbohm et all 120021) . it was then possible to com- 
pare these pulse-profiles and determine the shift in seconds be- 
tween any two profiles. Using these shifts, the starting period, 
and the number of elapsed pulses between the two profiles, we 
performed a linear fit to obtain a refined period and possibly a 
P with the refined period until no shift was detectable between 
the first profile at the start of the observation and the last profile 
at the end o f the observation (for a detailed explanation of this 
method, see lFritz eT al. 2006). With this technique, we were able 
to derive a refined averaged pulse period of 283.5320 + 0.0002 s. 
Our analysis, however, showed, that the pulse period exhibited 
fluctuations on timescales shorter than one orbital period with 
an amplitude of up to 0.003 s (for a detailed discussion of the 
evolution of the pulse period, see Staubert et al., in preparation). 
No pulse period "glitches" were observed, which is consistent 
with the short duration of the flares, during which no significant 
angular momentum transfer onto the neutron star is expected. 

The pulse-profile of Vela X-l is k nown for its remarkable 
energy dependence (see for example lLa Barbera et al.l 2003; 
iKrevkenbohm et al.ll2002l 119991: iRaubenheimerll 19901 and refer- 
ences therein) and, despite its complicated structure, also for its 
long-term stability (Rauben heimerlll990t IStaubert et~ai] [l980). 
While the complicated profile for energies below lOkeV con- 
sists of up to 5 peaks, the profile at energies >20 keV consists of 
two simple energy independent peaks. 

The pulse-profiles show n by lLa Barbera et all d2003l) ; 
Kretschmar et al. | (119971) . and IStaubert etalf Tl980). however, 
showed that the high-energy pulse profile of Vela X-l is also 
quite complex: around 30keV, the main peak shows a trian- 
gular shape, i.e. a linear rising flank followed by a sudden 
an d sharp fall while the secondary peak is sinusoidal (Fig. 3 
of iKretschmar et aTl 1 19971) . Pulse-profiles were obtained for all 
data including the flares, but without the eclipses (see Fig. [7]). 
The INTEGRAL data confirm these previous results and show 
that the secondary peak remains sinusoidal from 15keV up to 
100 keV, while the main peak evolves from a simple almost si- 
nusoidal shape around 20keV into a "triangular" shape with a 
smooth linearly rising flank followed by a steeply falling edge 
at 40keV and above. Apart from confirming earlier results, our 
INTEGRAL data therefore indicate that the high-energy pulse- 
profile is stable on timescales of decades and over significant 
luminosity ranges. 

We also derived pulse-profiles for flare 1 only. Despite the 
dramatic intensity increase, however, the number of pulses added 
was rather small. The signal-to-noise ratio diminished, and sig- 
nificant systematic artefacts became evident in the folded light 
curve due to the low number of pulses. In any case, no significant 
change in the pulse-profiles was evident within the uncertainties. 
We therefore do not discuss these profiles further here. 

3.3. The eclipse 

During this long observation of Vela X-l, two eclipses were 
observed. When examining Fig. [2] the ecl ipses appear to be 
far lon ger than the 1.69+0. 06 d reported by IWatson & Gri ffiths 
dl977l) . A detailed analysis of the "technical" circumstances, i.e. 
the INTEGRAL orbits and the perigee passages of the satel- 
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Fig. 7. Energy resolved pulse-profiles of Vela X-l using all data 
outside eclipse including giant flare 1 . The folding pulse period 
is 283.5320s (see Sect. I3.21 >. The profile has been adjusted such 
that the pulse minimum is at phase 0. The pulses are clearly de- 
tected also in the highest energy band up to lOOkeV. 



lite, revealed that the perigee passages of INTEGRAL occur 
directly before both observed eclipses of Vela X-l (see Fig. 13. 
Fortunately, there are a few science windows between the end 
of the perigee passage (i.e. the start of a new revolution) and the 
ingress into eclipse, such that it is possible to determine the start 
and end of the eclipse accurately. Since the optical companion of 
Vela X-l is a ma ssive B0.5Ib sup er giant with an extended vari- 
able atmosphere (ISato et al.ll 19861) . it is difficult to determine the 
ingress and egress times of Vela X-l in soft X-ray observations 
to high precision due to strong photoelectric absorption (see e.g. 




1 2 
2003 December 2 



Fig. 8. Close-up of the egress of Vela X-l from the first eclipse 
using a pulse averaged light curve of Vela X-l (time resolution 
of 283.5 s). The flux increases exponentially with time until a 
normal flux level is reached as indicated by the solid red line. 
The egress time is determined as the time when the count rate 
was no longer consistent with zero and is marked by a thick ver- 
tical line. 



IStuhlinger et al.l [20071 and references therein). At higher ener- 
gies, such as in the 20^-0 keV INTEGRAL data, photoelectric 
absorption poses less of a problem, although the determination 
of the egress or ingress time is only possible with an uncertainty 
of several 100 s (see Fig. [8), due to the changing extension of the 
atmosphere of the companion star and the general variability of 
the source. 

Using a barycentered pulse averaged light curve, we de- 
termined the ingress and egress times of both eclipses with 
an accuracy of about one pulse period (see Fig. [8). Using 
these times, we derived the duration of the eclipses in the 20- 
40keV energy ba nd to be 1.70 + 0.01 d, consistent with the 
1.69 d reported bv lWatson & Griffiths! d 1977b . Then we derived 
the center times of the two eclipses to MJD 52974.223 and 
MJD 52983. 195, respectively. Using these two times, we derived 
a new mid-eclipse reference time of MJD 52974.227+0.007. 
Note that the uncertainty in our measurement of the mid-eclipse 
reference time is far smaller than for earlier measurement s 
( Ivan der Kris & Bonnet-Bidaudll984tlWatson & Griffiths! 19771) . 
since the time resolution is significantly higher. To compare 
this mid-eclip s e time with the Tgo of the ephemeris given by 
iBildsten et all d 1997b . we need to convert the mid-eclipse time 
r ec i to the time of mean longitude T^q. Since the eccentricity 
e of Vela X- 1 is non-zero, this conversion is not straightfor- 
ward. The 790 is slightly offset from the center of the eclipse by 
AT = r ec i - Ti)Q (for a purely circular orbit, 790 and r ec i are iden- 
tical, and see also Fig. 0). This offset is given by dDeeter et al.l 
U9871) : 



Ar 



Port, L 1 (sin ;-/?)( 1 -I sin Q \ 

ccosw 1 + - r- + r (1) 

n \ 2 tan 2 i - - . 1 



2ySsin z / 



where j3 depends on the semi-major axis and the stellar radius: 



j6 = 



1 - 



(R/a) 2 
l-e 2 



(2) 



For the case of Vela X-l, p = 0.8 dDeeter et al.l Il987l) . 
Applying the corresponding values from Table |4] we obtain 
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Table 4. Ephemeris of Ivan Kerkwiik et ail £[995) and 
Bilds ten et al.l dl997|) . which have been used for obtaining 
the pulse period and the improved mid eclipse time r ec i, the 
corresponding Tgo, and the improved orbital period f or b. See 
text for discussion of the new Tgo. 



a) 



Vela X-l 



7*90 
Tgo 
T e d 
Tgo 

P orb 
P orb 
-^orb 

a sin i 



MJD 44278. 5466: 
MJD 48895. 2186: 
MJD 52974.227 : 
MJD 52974.001 : 
8.964416: 
8.964368 : 
8.964357 : 



: 0.0037 

:0.0012 

0.007 

:0.012 

:0.000049 d 
:0.000040 d 
:0.000029 d 



113.89 lt-sec 
i > 73 ° 

Ecc. e 0.0898 ±0.0012 

_oj 152.59 ±0.92 

References. (1) van Kerkw iik et alj <1995); (2) Bild sten et ail dl997h : 
(3) this work 



(1) 

(2) 
(3) 
(3) 
(1) 
(2) 
(3) 
(2) 
(1) 
(2) 
(3) 



a significant offset of 0.226+0.005 d (the uncertainty is due 
to the uncertainty in the inclination z). We therefore calculate 
MJD 52974.001+0.012 to be the new measurement of T 90 (in 
good agreement with the T 90 - 52974.008+0.008 obtained from 
pulse-timing analysis by Sta ubert et al., in prep.) . Comparing 
this Tgo with the ephemeris of lBildsten et al.l (1 19971) . we find that 
the time of mean longitude Tgo differs by about 0.005 d. This 
shift is unsurprising given that the uncertainty in the orbital pe- 
riod is of the order of 0.00004 d and the uncertainty in Tgo of 
iBildsten et all d 1997b was 0.0012d. Since then, Vela X-l has or- 
bited its companion 455 times such that the accumulated uncer- 
tainty amounts to a maximum shift of about 0.02 d, far more than 
the observed shift of 0.005 d. We therefore provide an updated 
Tgo and orbital period P 01D for the ephemeris in Table|4] but leave 
the other parameters of the ephemeris unchanged. 

3.4. Spectrum 

We combined all ISGRI data outside the eclipse to create a 
spectrum of high signal-to-noise ratio. Since a fully working 
physical model of the X-ray production mechanism in accret- 
ing X-ray pulsars does not exist due to the complexity of the 
problem (significant progress has, however, been made; see e.g. 
Becker & Wolff 2007), an empirical model has to be used, which 
consists of a power law modified by a high-energy cutoff at high 
energies. However, special care has to be taken not to introduce 
any artificial features by a "break" or a sudden onset of the cutoff 
as in t he case of the "High-energ y cutoff" in XSPEC (see for ex- 
ample Kreykenboh m et alj 19991 and references therein). This is 
true particularly if the source exhibits cyclotron resonant scatter- 
ing features (CRSFs). We therefore used a power law modified 
by the Fermi-Dirac cutoff (Tanaka 1986) to model the spectrum. 
This model describes well the continuum of Vela X-l in the 
ISGRI range, i.e. above 20keV, and was in fa ct used successfully 
in the past to fit the spectrum of Vela X- 1 ( Kreyken bohm et all 
119991) and other sources. To model the well known CRSF in 
the spectrum of Vela X-l at 50keV dKendziorra et all 1 19921) 
we us ed a Gaussian optical depth profile (GABS. ICoburn et all 
|2002|) . The full model for the ISGRI spectrum is then given by 



/«,„,(£) (x E T x 



1 



exp 



\ Efoid / 



+ 1 



Xexp(-T G ABS(£)) (3) 



A systematic error of 1 % was applied to account for the uncer- 
tainties in the response matrix of ISGRI. 
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Fig. 9. Spectrum of Vela X- 1 as obtained by INTFGi?AL-ISGRI 
a data and folded model, b residuals when using only a power 
law modified by the Fermi-Dirac cutoff, c residuals for the best- 
fit model with a cyclotron absorption line at 53.4^g keV. See 



Table [5] for all spectral parameters. 



Table 5. Fit results for the overall spectrum excluding the eclipse. 
When fixing r to one of the values found in the literature (fixed 1 
and fixed 2; iKrevkenbohm et alj|1999t lLa Barbera etaill2003l) . 
£ C ut and E F are better constrained, but the resulting -values 
are higher. All uncertainties here and elsewhere in the paper are 
on a 90% confidence level. The DOF are the degrees of freedom. 



Parameter 


free 


fixed 1 


fixed 2 


Exposure 


[ksec] 


508 


508 


508 


r 




1 6 +u 1 
lM -0.6 


1.8 fix 


2.0 fix 


Edit 
Efoid 


[keV] 
[keV] 


35 6 +7 5 
11.2 +0 * 


41.3+i-S 

10 9 + ^ 

-fl-3 


473+ Q6 
10 + " : * 
-fl.fi 


E c 


[keV] 


53 4 + ^ 


53 6 +8 - 6 
' -Q.6 


53 5 +ii -' 

-fl-3 


fX C 


[keV] 


76 ^ 
10* 

LM -OA 


1 •"-(). 1 


8 - 2+ 02 

1 l +!i: 2 

1 - 1 -o.i 


^ed(DOF) 




1.3 (20) 


1.4 (21) 


1.6 (21) 



The broad band continuum is described accurately by a 
power law with Fermi Dirac cutoff, but the spectral parame- 
ters were not well constrained (see Table [5]). This was expected, 
since no data below 20 ke V were available to enable us to deter- 
mine the spectral slope, the cutoff, and folding energy simultane- 
ously, especially since the cutoff energy is expected to be around 
~20keV, close to the lower end of the energy range. Fixing F 
to some (arbitrary) value from the literature allows to £ cut and 
Efoid to be constrained more accurately, however, depending on 
the choice of T also results in a rather bad fit (see Table [5j. In 
any case, the parameters of the CRSF do not depend strongly on 
the choice of F. 
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After fitting the broadband continuum, highly significant fea- 
tures remain, which are due to the cyclotron line at ~50keV 
(see Fig. [9J3). After the inclusion of a Gaussian absorption line 
at 53.4*91 keV, the resulting fit is acceptable = 1.3 with 
20 degrees of freedom) and no significant features remain (see 
Fig. |9]and Table[5]l. Note that the shallow line like residuals be- 
low ~30keV in Fig. [9j} are not due to the disputed CRSF at 
25keV but merely a consequence of the 50keV line. The ques- 
tion of whether the 25 ke V line exists or not cannot be answered 
by this observation because no data below 20keV are available, 
which would be crucial to determine the continuum and to de- 
tect a CRSF at 25 keV0. The resulting best fit is shown in Fig. [9^ 
and the residuals in Fig [9};. Note that after the inclusion of the 
Gaussian absorption line at 53.4keV, some residuals remain (see 
Fig- lit) at 80keV. In this energy range, some strong background 
lines are present i n the ISGRI background (tungsten and lead, 
iTerrier etaL 2003) that might be responsible for these residu- 
als. Ho wever, it is well known that CRSFs have a complex line 
shape dSchonherr et alj l2007t lArava-Gochez & Hardin J 120001: 
lArava & Harding! 119991) and are therefore not well modeled 
by simple Gaussian or Lorentzian functions (s ee for exam- 
ple iKrevkenbohm et al l l2005t iHeindl et al.lll999l) . The residu- 
als present in Fig. [9J; could therefore also be due to an insuffi- 
cient description of the CRSF by a Gaussian absorption line and 
therefore incorrect modeling of the underlying continuum. To 
improve this unsatisfactory situation and to derive real physical 
parameters from the line shapes, effort s are made to create li ne 
shapes using Monte Carlo simulations dSchonherr et al. 2007). 

A more detailed analysis of the spectrum using various spec- 
tral models, studying the evolution of the spectral parameters 
with time, and, in particular, using phase resolved spectroscopy, 
and a comparison with previous work is beyond the scope of this 
paper and will be discussed in a forthcoming publication. 



3.5. Spectral evolution during the flares 

After studying the averaged spectrum, we now consider the spec- 
tral evolution of the source during the flares, especially giant 
flare 1 . Figure|2]3 shows the hardness ratio over the flare, defined 
to be 

HR=^ (4) 
H + S K ' 

where H is the count rate in the hard band (40-60 keV) and S 
the count rate in the soft band (20-30 keV). The hardness ratio 
remained constant throughout most of the observation, i.e. no 
correlation with orbital phase was evident. During the flares, the 
behavior of the hardness ratio, however, was remarkable. Shortly 
before the onset of giant flare 1 as well as during the flare, a clear 
deviation in the hardness ratio from its overall average value of 
-0.74 was evident; the hardness ratio declined suddenly to -0.82 
and later during the flare even to -0.85 (see Fig. [2J3). A similar 
behavior was observed for flare 4: with the sudden onset of the 
flare, the hardness ratio dropped from a pre-flare value of ap- 
proximately -0.72 to -0.84, the same level as in giant flare 1, 
although flare 4 was far shorter and reached only a third of the 
peak flux of flare 1. This result means that the source became 
significantly softer during these flares. Strikingly, however, this 
softening did not always occur: during flares 2 and 5, no sig- 
nificant change in the hardness ratio was apparent. In fact, the 
hardness ratio during these flares remained constant at the av- 
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3 Using INTEGRAL data from a 2 Msec long observat ion in 
INTEGRAL A03 including JEM-X data, dSchanne et aT]|2006t) , how- 
ever, confirmed the existence of the 25 keV CRSF. 



Fig. 10. Hardness Intensity Diagram of Vela X-l; data from the 
eclipses have been excluded. The time resolution is one science 
window (typically 1800 s). The hardness ratio is defined as given 
in Eq.[4] The data points from the flares (see Fig.|2]and Tabled 
are indicated by individual symbols. Note that flares 1 and 4 both 
show significant softening opposite to flares 2 and 5, which show 
no softening (see Sect. 13.5b . No data are available for flare 3 (see 
Sect. ED- 



erage value. This result is remarkable, since with a duration of 
~9 h, flare 5 lasted far longer than the other flares and its peak 
flux reached >50% of the peak flux of giant flare 1 ; it was there- 
fore far longer and brighter than flare 4, but it is associated with 
no spectral softening. 

Although the source was bright during flares 1 and 5, it was 
not possible to derive meaningful spectral fits for spectra from 
both flares since the net exposure time was still short. We refrain 
from discussing the spectra of the flares in more detail since the 
spectral parameters could not be reliably constrained. 

"Hardness intensity diagrams" (HIDs) are a useful way to 
study the spectral evolution of a source. A HID for the full obser- 
vation of Vela X-l, excluding the eclipses, is shown in Fig. [10] 
As expected, the HID of Vela X-l does not show the q shape 
observed for black holes dFender & Bellonill2004l) . In Vela X-l, 
most of the data points are clustered around the average values 
of intensity and hardness ratio. The only exception are the data 
points from giant flare 1 and the other flares, which are above the 
general cluster of data points due to their - by definition - higher 
intensity. The data points of flare 1 and flare 4 are shifted, due 
to the spectral softening during the flare, as is already evident 
from the hardness ratio (see Fig. [2)). Note that the softening did 
not evolve during the flares, but the flares are softer than the av- 
erage spectrum from the very onset. On the other hand, flares 2 
and 5, although of comparable intensity, did not show any soften- 
ing. Instead, the data points of the flares form two distributions 
(flares 1 and 4 versus flares 2 and 5) that are disjunct (see Fig.fTOt. 
On the whole, the HID shows a slight hardening towards lower 
count rates. 

Another method of searching for spectral evolution are color- 
color diag rams (CCDs); e .g. V 0332+53 shows an interesting 
evolution (Rei g et al.ll2006l) . We derived a CCD for Vela X-l but 
no interesting behavior was found. 
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4. Discussion 

Vela X- 1 has been known for a long time to have a highly time- 
variable light curve and to show intensity variations of up to a 
multiple or a fraction of the original intensity on time scales of 
days, hours, or even seconds. 



4. 1 . The flares 

Although Vela X-l has exhibited extensive flaring activity in the 
past (see iKrevkenbohm et al. ll999tlHa berl 1994; Lapsho v et alj 
I l992t lHaberl & Whitelll990t lNagase et alJ I 1983. among others), 
however, giant flares (as flares 1 and 3) had not been seen be- 
fore. The first question is therefore whether these flares are a 
rare phenomenon or just could not be detected. The problem 
is that all these flares are rather short - even giant flare 1 and 
long flare 5 have a duration of a few hours. Unless the source 
is constantly observed as in our observation, typical sky moni- 
tors such as the All Sky Monitor (ASM) onboard RXTE fail to 
detect these flares. The RXTE-ASM monitors the entire observ- 
able sky regularly, typically observing Vela X-l up to 10 times a 
day, but often only two or three 90 s long dwells per day are avail- 
ably This monitoring (which is furthermore irregular) does not 
allow us to de tect such short-lived events as these giant flares. 
As shown by IStaubert et all (120041) . the RXTE-ASM failed to 
detect the giant flares 1 and 3 when only taking 90 s dwells into 
account. During flares 2 and 5, it showed only minor increases, 
but given the typical uncertainty in the ASM data and the over- 
all scatter, no unusual behavior of Vela X-l could be detected. 
We note that flare 1 was present when relaxing the assumption 
of using only 90s ASM dwells to dwells with >80s exposure; 
however, these dwells often have larger ;^ 2 ed -values than recom- 
mended. Furthermore, detecting short-term variability with the 
ASM is problematic, since the ASM sometimes also erroneously 
detects significant flux (even at the level of the flares) while 
Vela X-l is in eclipse. We therefore conclude that monitoring 
flares on a single dwell level with the ASM is unreliable and 
uninterrupted observations are required. This finding also read- 
ily explains why these giant flares have to date been unnoticed, 
since long uninterrupted observations spanning at least one full 
binary orbit of Vela X-l are only rarely performed. 

The analysis of the hardness ratio of the flares in Sect. 13.51 
shows that there seem to be at least two different types of flares 
(see also Fig.lTOli: the first type (flares 1 and 4) showed dramatic 
increases in the count rate and the onset of the flare was very 
sudden. These flares showed significant spectral softening dur- 
ing the flare. They could appear at any time, i.e. were unrelated 
to the previous evolution in a similar way to flare 1, which is 
superimposed on a downward trend. Although flare 4 was short, 
the change in the spectral hardness ratio during the flare was sig- 
nificant (see Fig. |2j)). The second type was more similar to a 
high intensity state than an actual flare. Flares of this type are 
longer than the first type. The hardness ratio does not change 
during these flares, indicating that the spectrum does not change, 
and the source simply becomes brighter. For the short flare 3, 
no hardness ratio could be obtained since INTEGRAL was in 
engineering mode during that time; but given that it is short and 
features a very dramatic rise, it is likely that it belongs to the first 
type. 

Although the hardness ratio in Fig.[2j3 is already conclusive, 
a hardness ratio comparing the high-energy spectrum (above 
20keV) with the low-energy spectrum (below lOkeV) would 
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be even more interesting since changes due to photoelectric ab- 
sorption are only visible at low-energies. However, no JEM-X 
data are available. It is therefore impossible to determine how 
bright the flares were at lower energies with our current data set. 
Given the flux level of more than 5 Crab in the peak of flare 1 
and the significant softening observed (see Sect. 13.5b . we can as- 
sume that the source was also extremely bright at lower energies. 
However, high photoelectric absorption could have dampened 
the brightness again in the classical X-ray band. The evolution 
of A^h during a flare could possibly help to differentiate between 
flare types and their underlying mechanisms, which is impossi- 
ble with our current data set. 

The mechanism behind these different types of flares, how- 
ever, is not understood. Simulating asymmetric adiabatic ac- 
cretion flows onto a neutron star in the wind of OB stars, 
Taam & F rvxelll d 19891) demonstrated that a temporary accretion 
disk may form in systems such as Vela X-l. The formation of 
this temporary accretion disk is the result of an interaction be- 
tween the incident flow and shocks in the wake region and is a 
general property of a binary system consisting of a neutron star 
and an OB companion. Another consequence is a reversal of the 
accretion flow. Associated with the flow reversal is the destruc- 
tion of the accretion disk, resulting in a significantly increased 
accretion rate. During this short phase, the material stored in 
the temporary accretio n disk is accreted onto the neutron star. 
iTaam & Frvxelll (119891) showed that these flow reversals occur 
on timescales of several hours, and predicted flares that would 
last from 15 to 60 minutes, in agreement with the observation 
of flares 3 and 4 and also with the 1 h long flare observed by 
iKrevkenbohm et aTl d 1999b . Furthermore, the overall flaring re- 
cu rrence timescale durin g our observation agrees with that given 
bv lTaam & Frvxelll (Il989b . 

Later hydrodynamical studies found that wind-accretion 
onto the neutron star is already a highly instable process in itself. 
The accretion wake following the neutron star contains dense fil- 
aments of compressed gas with density variations of a factor of 
100 compared with the undisturbed wind. When accreted, these 
density fluctuati ons produce abrupt changes in the X-ray lumi- 
nosity (see e.g. iBlondin et al.l 11990b . Even more important for 
wind-accretion, however, is the velocity structure of the wind, 
since 

Lx oc 4 (5) 

where p is the dens ity and v the velocity of the wind 
dBondi & Hovld 1 1 9441) . Hydrodynamical simulations have 
shown that not only the density but also the velocity of the stellar 
wind changes d ramatically with time including sharp drops and 
spike s (see e.g. lRunacres & Owockill2002t iRunacres & Owockil 
2005). Furthermore, the shock trailing the neutron star oscillates 
with brief periods of disk formation, forcing the accretion flow 
to ch ange its pattern, generating the so-called "flip-flop instabil- 
ity" dMatsuda etal.lll99ll Il987h . This instability then produces 
disk-l ike rotational inflows t hat change their direction repeatedly 
(see [Benensoh n et al.lll997l and references therein), and is an 
intrinsic phenomenon that occurs whenever a gas stream flows 
past a neutron star or black hole, no dens ity or velocity gra- 
dient being necessary dMatsuda et al.lll99ll) . The timescale for 
this flip-flop behavior is calculated to be of the order of 45 min 
(Benenso hn et al.l 1 1997b. which wou ld agree with the flare ob- 
served by Rrevkenbohm et al.1 d 19991) and flares 3 and 4 in Fig. [2] 
On the other hand, the flip-flop behavior does not explain flare 1, 
which is superimposed on a general downward trend, because 
it is far longer (several hours), nor flare 5 which is about 12 h 
long, since the calculated accretion rates vary on an even shorter 
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times cale of ~100s or less (see e.g. Fig. 3 of iBenensohn et al] 
fl997h . 

While the above scenario explains well the short flares seen 
during our observation (flares 2, 3, and 4) and the overall rapid 
variability of Vela X-l, the flip-flop instability fails to explain 
the intense long flares presented in this paper (flares 1 and 5) 
that last several hours; another mechanism must therefore be at 
work that can alter M by a factor of up t o 10, not only for ~100 s, 
but for many hours. iKaper et alj d 19931) demonstrated that, due 
to the dumpiness of the shocked wind, the local density varies 
by a factor o f 100, which can expl a in the flaring X -ray luminos- 
ity (see also lOskinova et alj|2007h . iLevder et al.1 d2007l) argued 
that dense clumps trapped in an otherwise thin and more ho- 
mogeneous wind might be responsible for long flares, when the 
clumps are being accreted. Since giant flare 1 lasts several hours, 
however, it is unlikely that a single blob could feed the accretion 
for such a long time, given that the 6 h and 12 h durations of giant 
flare 1 and of flare 5 correspond to a considerable fr action of the 
neutron star orbit. Instead, t he clumpy OB star wind ( Kap er et ail 
Il993t iLucv & White! [19801) is probably viscously smeared out 
in a (small) accretion disk of the system. This filled accretion 
disk can then feed the neutron star with a significantly higher 
M than usual over several hours. A change of P, however, can- 
not be observed, since the transferred angular momentum during 
such a flaring episode is far too small (see also Sect. I3.21 >. When 
Vela X-l is less active, the OB star wind is probably less struc- 
tured. In summary, we conclude that the observed long flares are 
due to a strongly structured OB star wind. 

Since it is difficult to detect or monitor the flares unless 
during a pointed observation, it is not possible to determine 
how frequent such events are. Apart from the flares in our ob- 
servation, reports of flares — although far less inten s e - are 
common in the literature dKrevkenbohm et al.l 119991 : [Haberl 
Il994t lHaberl & Wriltel 1 19901 among others). In another long 
(~2 Msec) observation of the Vela X-l region in 2005 November 
by INTEGRAL, Vela X-l again exhibited very intense flares 
dKrevkenbohm et alj 120061 ISchanne et al l 12006). We therefore 
conclude that bright flares are quite common in Vela X-l. 
However, we caution that Vela X-l does not always show high 
activity. During a long INTEGRAL observation in Summer 
2003, Vela X-l was in a quiet phase with only little flaring ac- 
tivity detected. 



4.2. Short term variability 

Although the flip-flop instability discussed in the previous sec- 
tion fails to explain the large flares exhibited by Vela X-l, 
it could explain its short-term variability: Vela X-l exhibits 
strong pulse-to-pulse variations and other short-term fluctua- 
tions superimposed on the general variability, while at the same 
time the long-term pu l se-profile remains constant o ver decades 
dRaub enheimerl Il990t IRrevkenb ohm et aT] 120021) . iTaam et ail 
d!991l) developed several scenarios and the resulting behavior 
of the mass accretion rate. These authors obtain a time variation 
in M that exhibits short flare-like events on time scales of ~ 100 s, 
i.e. less than one pulse period of Vela X-l, for an accretion rate 
of 0.7Mhl, where Mhl is the H oyle-Lyttleton mass capture rate 
in terms of the Eddington rate dHovle & Lvttleton|[l94ll 119391) . 
and a wind velocity of about 1000 km s ; such behavior is typ- 
ical for O- and B-type stars and precis ely that required to ex - 
plain the observed short-term variability. Watanabe et al. (2006) 
studied the stellar wind of the Vela X-l system and were able to 
explain the observed line inte nsities by using a wind model de- 
veloped bv lCastor et alj d 19751) with a terminal wind velocity of 



1100 k ms which agrees with the calculations of ITaam et al.l 
(fl99l . Since this instability and the transfer o f angular mo- 
mentu m appear to be intrinsic to wind-accretion (Matsuda et al. 
1991), the accompanying torque reversals can therefore also be 
expected to exhibit randomly varying short spin-up / spin- down 
episodes, or a "flip-flop" behavior. As lBovnton et alj £l984) sug- 
gested, this instability might therefore be identified with the 
random-walk behavior of the pulse period, observed in many 
wind-accreting sources including Vela X-l. 



4.3. The off states 

In a simil ar way to t he fla ri ng activity, the off stat es re- 



ported bv llnoue et alj d!984l) . iKretschmar et all dl999h . and 
Kreykenb ohm et aiT i 1999l) r~where the source was not simply 
weaker but below the detection limit of the respective instru- 
ments and no pulsations were observed, are remar kable. After 
the off state observed by IRrevkenbohm et aT] d 19991) . the source 
resumed its normal , pulsating behavior wit hout any transition 
phase (see Fig. 4 of lKjev kenboh met al.lll999h : also no unusual 
behavior of t he sou rce was observed after the end of the off state. 
Ilnoue et al] d 19841) reported that these off states occur without 
any prior indication. In our observation, we observed several 
short off states (see Table |3j, which also occured without any 
prior indication and without a transition phase such as a slow 
decay. Furthermore, as shown in Fig. [5] off state 1 occurred dur- 
ing a phase in which the source had an average intensity level 
of about 250 mCrab. The count rate decreased dramatically to 
zero for ~850 sec (corre sponding to three pulse peri ods) with no 
pulsations being visible. Kreykenb ohm et all d 1999b reported ob- 
serving Vela X-l in an off state of 550 sec (corresponding to 2 
pulse periods) at the start of an observation, although the total 
duration of this off state is unknown, and throughout the obser- 
vation, no p ulsations were detected. It i s possible that the events 
observed bv lKrevkenbohm et al.l dl999) and in Fig.|6]are similar, 
i.e. the observation of the former started just in the middle of a 
short off state similar to that in Fig. [6] 

The reasons for these off states and the sudden reappearance 
of pulsations are not understood, but Vela X-l appears to expe- 
rience them on a regular (but non-periodic) basis (see Fig. [6] 
and Table [3). Although several explanations of thes e phenom- 
ena h ave been put forth (including transiting planets, Hav akawal 
11984 . none can explain the observed off states, since all of these 
ideas (for example, variations in the mass-loss rate of the optical 
companion) req uire a significantly lon ger timescale. "Blobs" in 
the stellar wind dFeldmeier et al.lll997l) . for example, would not 
only require extremely high optical depths to block the X-rays 
and gamma-rays completely, but would also need very sharp bor- 
ders to explain the observed sudden turn-off and turn-on of the 
source (sometimes within a single time bin of 20 s, as shown 
in Fig. |6j. We therefore consider it highly unlikely that typical 
clumps in the stellar wind are responsible for the off states shown 
in Fig. [6]and other explanations must be considered. 

The wind of early-type super giants like HD 77581 is known 
to be inhomogeneous and clumpy (see discussion above and both 
IWalter & Zurita-Herasll2007t iBlondin et al.lll990h . Typical mod- 
els show that the density in the stell ar wind of super giants can 
vary by several orders of magnitude dRunacres & Owockil2 005). 
These density variations should not only correspond to the pres- 
ence of clumps - i.e. regions of strongly increased density - but 
also holes - regions of strongly reduced density. In these holes, 
the density is l ower than the average den sity of the wind by 
a factor of 10 3 dRunacres & Owoc ki 2005). If the neutron star 
enters these holes, M would then also decrease by a factor of 
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~ 10 3 and the X-ray luminosity would be reduced accordingly. 
Furthermore, the density fluct uations predicted by these m odels 
occur suddenly (see Fig. 1 in iRunacres & Owockill2005l) . as is 
the onset of the off states (see Fig. [2j. On the other hand, these 
models predict that the density always varies; it could therefore 
be expected that the off states should be rather common, which 
does not seem to be case. Therefore, some additional mechanism 
must be at work that is triggered only rarely. 

Another me chanism that can explain the se off states is the 
propeller effect (llllarionov & Sunvaevlll975l) . In short, the pro- 
peller effect inhibits the accretion of material onto the compact 
object when the Alfven radius (where the ram pressure of the in- 
falling gas and the magnetic pressure are equal) is larger tha n the 
co-rotation radius dPringle & Reeslll972tlLamb et al.lll973"t) . The 
Alfven radius, however, is not constant: it depends on the amount 
of infalling material, M. If M drops, the Alfven radius will in- 
crease and once larger than the co-rotation radius, accretion is 
no longer possible and the X-ra y source ba sically switches off 
and no pulsations are observable . ICuil dl 997 ) observed this effect 
in GX 1 +4: in very low luminosity states, no pulsations were ob- 
servable while in high luminosity states, the source was strongly 
pulsating (however, other explanations are also possible to ex- 
plain the absence of pulsations in GX 1 +4; see Ferrigno et al. 
I2007L and references therein). The magnetic field strength for 
which the system enters the propeller regime is then given by 
(ICuil l 1997b 



B = Cx 



P \ 7/6 



(n) 



10~ 9 erg crrr 2 s _1 \lkpc/ \1.4M, 



M 



1/3 



(6) 



where the constant C = 4.8 10 10 G, P is the spin period of 
the neutron star, F x is the bolometric X-ray flux, d is the dis- 
tance, and M is the mass of the neutron star. Since the strength 
of the magnetic field of Vela X-l is known f rom the observation 
of the cyclotron resonant scattering features ( Kreykenbohm et all 
2002), Eq.0 can be used to obtain directly the critical flux limit 
for the onset of the propeller effect: 



' X.Propeller 



4.3 x 10 7 erg cm 2 s 1 



: (i0 12 g) (is) 



p r 773 



1 kpc 



M 



1.4 M f 



-2/3 



Applying the corresponding values for Vela X-l, i.e. 
5 = 2.6 10 12 G, P = 283.5 s, d = 2 kpc, and M=1.9M 
(see Sect.[TJ, we obtain a critical flux of 



-Fx.Propeller, Vela X-l ~ 1-1 x 10 



-12 -2 - 

erg cm s 



Compared with the typical bolometric flux of several times 
10~ 9 erg cm -2 s -1 , this critical flux is lower by about three orders 
of magnitude. This observed bolometric flux corresponds to an 
intrinsic luminosity of ~ 6 x 10 32 erg s . We can therefore safely 
conclude that the propeller effect alone cannot force the source 
to switch off when Vela X-l is in normal accretion mode. 

A possible scenario would, however, be to combine both - 
individually unsuccessful - mechanisms: the neutron star was 
evidently in a region in the stellar wind of strongly variable den- 
sity, as demonstrated by the significant overall variability of the 
source and the presence of many flares. The models for the stel- 
lar wind predict density vari ations of up to a factor of 10 3 ~ 5 
dWalter & Zurita-He ras 2007). If the neutron star enters this re- 
gion with very low density, its luminosity will drop accordingly 
because the X-ray luminosity depends linearly on M, and pro- 
duce a luminosity of the order of less than ~ 10 33 ergs s _1 . As 



shown above at such low luminosity levels, Vela X-l enters the 
propeller regime. This explains why no residual pulsations are 
observed during the off states. 

Off state 5, which is considerably longer than the other off 
st ates and does not show a sudden end similar to the dip observed 
bylRretschmar et al.l d!999t) . exhibits, however, a different behav- 
ior. Unlike the off states shown in Fig. [6] these dips do not show 
a very sudden onset or end; instead their onset (and end) is rather 
smooth and pulsations could also be o bserved for some time af - 
ter the onset of the dip (see Fig. 2 in iKretschmar et all 1 19991) . 
These authors also observed a dramatic increase in the photo- 
electric absorption (more thanlO 24 cm -2 , Compton thick). These 
dips therefore belong to a second class of dips whi ch are caused 
by sig nificantly increased photoelectric absorption (ICharles et al.l 
1978) due to an optically thick cloud in the wind of the stel- 
lar companion passing through the line-of-sight. Such clouds 
are know n to be quite common in super-giant s ystems (see for 
exam ple iNagase et al.|[l98 6) and also Vela X-l (Watanab e et al.l 
2006). The existence of these clouds (also referred to as blobs 
or clump s) in the wind of an OB stellar companion was pro- 
posed by Luc y & White! (119801) . who suggested that the winds 
of OB stars themselves are not homogeneous, but break into 
a population of blobs that are radiatively driven through an 
ambient gas. This model was later modified to include radia- 
tively driven shocks in the stellar wind. Since then, models of 
OB star winds have u sually included blobs and shocks (see e.g. 
IFeldmeier et al.l 19971 and referenc es therein). During the dip ob- 
served by Kretschmar et al. (1999), the column density, TVh, was 
of th e order of 10 24 cm' 2 , similar to that measured by, for exam- 
ple, |LeyderetaTJ (|2007|) for the density of clumps in the wind of 
HD 74194/IGRJ08408-4503. 

We therefore conclude that off states, characterized by a sud- 
den onset and end, are usually short, and could be caused by 
a sudden drop in M that allows Vela X-l to enter the propeller 
regime. Intensity dips, however, are significantly longer, show a 
smooth transition, and spectra taken during the dip provide mea- 
surements of photoelectric absorption of more than 10 24 cm" 2 . 
These dips are readily explained by a dense blob in the wind 
passing through the line of sight. 



( ? ) 4.4. Connection with SFXTs 



The similarity between the flares and off states in Vela X-l 
and the behavior o f Supergiant Fast X-ray Transients (SFXTs, 
SgueraetalJ[2005) is intriguing. SFXTs are high mass X-ray 
binaries that show very brief outbursts on timescales of hours 
or even only tens of minutes , and then remain un detectable 
for months between outbursts (Negue ruela et al.ll2008l) . The gi- 
ant flares of Vela X-l reported in this work are similar to 
these outbursts, which are assumed to be due to the accretion 
of a dense blob of material embedded in a thin stellar wind 
dWalter & Zurita-He ras 2007): this would then imply that Lx is 
a direct tracer of the density of such blobs in the stellar wind. 
The same holds true for Vela X-l, for which giant flares such as 
flare 1 are also probably due to the accr etion of a den s e blob 
in the stellar wind, as discussed above; Ferri gno"et alj (2008) 
reached a similar conclusion in explaining the flares observed 
in IE 1 145.1-6141. However, we also conclude that short flares 
and the general variability can be explained well by the flip-flop 
instabi lity, and no additional cloud s in the stellar wind are nec- 
essarv. lGrebenev & Sunvaevl d2007l) considered, why SFXTs are 
unobservable when they are not in outburst. According to these 
authors, SFXTs should be bright persistent objects since the neu- 
tron star is deeply embedded in the stellar wind. Similar to our 
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explan ation of the off states of Vela X-l, iGrebenev & Sunvaevl 
(2007) invoke the propeller effect to explain the absence of de- 
tectable X-rays from the SFXTs when the sources are in quies- 
cence. Unlike Vela X-l, however, SFXTs are usually in the off 
state, while Vela X-l is in a normal accretion mode. The rea- 
son why SFXTs are switched off by default could either be a 
significantly thinner stellar wind (resulting in a lower M) or a 
significantly stronger magnetic field, which would inhibit the ac- 
cretion even at the densities encountered typically in the stellar 
wind. We therefore conclude that SFXTs and Vela X-l are very 
similar systems, except that SFXTs are normally in the propeller 
regime, while Vela X-l is normally in the accreting regime. Both, 
however, can switch sides, i.e. SFXTs can go into outburst, while 
Vela X-l can enter the propeller regime. 

4.5. QPOs 



Sect. 14. 3K as occured in off state 1 (see Fig. [5]). This scenario is 
interesting because it helps to understand the other off states in 
the following part of the light curve (see Fig. [5j, and to explain 
the density variations discussed in Sect. 14.31 

The quasi-periodic modulations in the source intensity could 
therefore easily be identified with similar structures in the wind. 
It should then be possible to draw conclusions from the period 
and observed intensity variations, about the applicable models 
and parameters and the underlying physics. However, these mod- 
els, as currently present in the literature, are general and appli- 
cable numbers cannot readily be ascertained. In principle, how- 
ever, as soon as these numbers can be obtained, the systematic 
and long-term probing of the stellar wind structure by the neu- 
tron star in Vela X-l could provide physical constraints of the 
various wind models. 



Several accreting X-ray pulsars show one or more quasi- 5. Summary 



perio dic oscillations (QPOs) in addition to the normal pulse pe- 
riod ( Shirakawa & Lai 2002). Accreting X-ray pulsars that ex- 
hibit long period QPOs include 4U 01 15+63 ( Pqpq ~ 500 s 
Heindl et alJI 19991) and V 03 32+53 (P QP0 ~ 20 s IMowlavi et all 



2006; Takeshima et al. 1994) among others. The reason for the 



existence of QPOs is not always clear, although the 0.05 Hz QPO 
in V 0332+53 i s thought to be due to inhomogeneities in the 
accretion disk (iMowlavi et ail 120061) . while a second QPO at 
0.22 Hz (lOu et al.ll2005F is almost (but not entirely) coincident 
with the pulse period of the system. It can therefore be assumed 
to originate in the X-ray production region, and maybe also be 
linked to the strong pulse-to-pulse variations in that system. In 
4U 1907+09, another wind-accretin g system, a transient QPO 
was observed with a period of ~ 1 8 s Jin 'tZand et alJ[T9 98). Due 
to the close similarity of this system with Vela X-l, we could 
expect to observe a QPO with a period between 10 s and 40 s. 
As discussed above, no evidence for a QPO with a period below 
140 s was found. However, a QPO with a long periodicity in the 
range of several thousand seconds was observed instead. 

The quasi-periodic behavior observed in Vela X-l (see 
Fig-HJ) is intriguing, albeit quite short lived. Such a periodicity - 
if real - in the luminosity of the source would indicate that the ac- 
cretion rate onto the neutron star also varies periodically, which 
implies the presence of a quasi-periodic structure in the stellar 
wind of the companion star, either in density or in velocity or 
both. Such periodic structures in the wind are certainly not due 
to the flip-flop instability discussed above, because the instabili- 
ties produce a chaotic non-pe riodic modu lation in the accretion 
rate (see e.g. Figs. 4 and 6 of iTaam et alJll99ll) . These periodic 
structures could be driven by instabilities in the atmospheres of 
very luminous stars, which are likely to reach distances of a few 
stellar radii. Various mechanisms for the production of these in- 
stabilities are discussed in the literature, amo ng them radiatio n- 
pressure-driven hydrod ynamical instabilities dShaviv 200TalbT), 
rotational instabilities (Fullerton et al. I ll997h . non-radial pulsa- 
tions dOwocki & Cranmerl |2002|) . and surface magnetic fields 
(IUd-Doulall2002l) . 

One of the phenomen a predicted by some models (e.g. 
lUd-Doula & Owockill2002h are "ray" like structures in the stel- 
lar wind due to the magnetic field of the OB star. It can be imag- 
ined that when the neutron star passes through these more or 
less periodic ray structures in the wind, the varying mass accre- 
tion rate M could then lead to the quasi-periodic behavior shown 
in Fig. [5] When the neutron star enters a region where the stellar 
wind has a significantly lower density, the accretion rate M could 
drop to a level at which Vela X-l enters the propeller regime (see 



We have presented the analysis of a long, continuous 
INTEGRAL observation of Vela X-l spanning about 1.7 binary 
orbits. In detail, our results are: 

- Vela X-l was found to be in a highly active state; 

- the observation of several flares were recorded, two of them 
giant flares with a peak brightness of more than 5 Crab; 

- a transient QPO with a period of ~6820s was observed, 
which we attributed to inhomogeneities in the stellar wind; 

- several off states were identified during which the source be- 
comes undetectable by INTEGRAL and no pulsations were 
visible. The onset of the off states was very sudden, in addi- 
tion to their end; 

- the pulse period was determined to be P = 283.5320 ± 
0.0004 s with no evidence for a spin-up or spin-down during 
the entire observation; 

- non-sinusoidal high-energy pulse-profiles were obtained up 
to lOOkeV; 

- the spectrum exhibited the well known CRSF at 53.4 keV, but 
due to a lack of low-energy data, the elusive line at ~25 keV 
could not be observed; 

- Vela X-l exhibited two types of flares: rapidly rising flares, 
which were correlated with a spectral change, and high inten- 
sity states, which were longer and during which the spectrum 
remained unchanged; 

- the short flares can be explained well by the flip-flop instabil- 
ity: the predicted timescales of 15 to 60 minutes are in perfect 
agreement with the observed duration of the short flares; 

- two different types of off states exist, which are probably ei- 
ther caused by a dense blob blocking the line-of-sight, or the 
onset of the propeller effect due to a drop in M; 

- the similarity between Vela X-l and SFXTs is striking: giant 
flares (in SFXTs: outbursts) are probably caused by the ac- 
cretion of a dense blob from the stellar wind, while off states 
(in SFXTs: quiescence) are likely to be caused by the onset 
of the propeller effect. 
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Fig.A.l. Dependence of the Crab count rate as measured by 
HJight on the off-axis angle. The average count rate of the Crab 
when the source is in the fully coded field-of-view (4!5) as ob- 
tained by the standard IBIS pipeline is shown by the dashed line. 
Note that count rates obtained by HJight are within 5% of this 
nominal value even up to high off-axis angles. HJight, however, 
almost always slightly underestimates the count rate compared 
to the standard pipeline. 
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Appendix A: Performance of HJight 

For our analysis of the temporal properties of Vela X-l, we require a bin size of 
20 s. The standard IBIS analysis pipeline, however, can only create light curves 
with a minimum bin time of at least 60 s and typically a bin time of 100 s is used. 
We therefore used the tool HJight (part of OSA 7.0) to derive light curves, since 
it is not only efficient and convenient to use but has the capability to derive light 
curves with a time resolution as high as 0.1 s. Since the IBIS cookbook cautions 
that HJight should only be used to analyze the temporal behavior within a given 
science window, we carefully verified the performance of HJight by using all 
publicly available Crab data up to revolution 464, i.e. about 1500 science win- 
dows, corresponding to an exposure time of about 3 Msec. First, we operated the 
standard IBIS pipeline up to the imaging level to derive background and off-axis 
correction maps to correct for the absorption due to the NOMEX effect (for expla- 
nation, see the IBIS analysis manual, and for the implications on the calibration, 
see lLubihski et alJ2 004); no jumps in the count rate with changing off-axis angle 
were therefore observed. We then used HJight to obtain light curves. Using the 
same time resolution of 20 s as for the Vela X- 1 analysis, we derived an average 
count rate for the Crab of 172.5 counts s between 20keV and 40keV when the 
source was in the fully coded field-of-view (max. off-axis angle of 4?5), which 
is very close to the 176.8 counts s obtained by the ibis pipeline for an entire 
science window. Furthermore, we found that the deviation in the count rates ob- 
tained with HJight is within five percent of the nominal value up to off-axis 
angles of more than 14°, which is of the order of the expected uncertainty (see 
Fig. IA.lt . In any case, HJight systematically underestimates the count rate (see 
Fig. IA.lt . and the derived flux values can be considered lower limits. In sum- 
mary, we found that HJight produces usable count rates up to an off-axis angle 
of 14°, we therefore conclude that the usage of the HJight tool for our purposes 
is well justified. 
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Fig. 2. Variability of Vela X-l for the complete Vela region observation from Revolution 137 to 141. a ISGRI 20^-0 keV light curve 
(time resolution 1 SCW, i.e. ~1800s) and b 20-30keV vs. 40-60 keV hardness ratio as defined by Eq. H] (rebinned by a factor of 
~3 with respect to the light curve). Labels indicate the revolution number. Short vertical lines below the X-ray light curve show 
INTEGRAL'S perigee passages, during which the instruments are switched off. The long dashed vertical lines show the ingress and 
egress times. The dotted vertical line indicates the derived eclipse cente r (see also Table|4|, while the dash-dotted line indicates the 
time of mean longitude Tgp based on th e ephemeris from lNagasd (Il989l) . Note that the offset of the newly derived (see Tabled 
in comparison to that of Nagase (1989) is too small to be visible in this figure. See text for further discussion. 



